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Imaging method e.g. for neuronal connectivity reconstruction, determines 
time variation in matter concentration for image elements of successive 
images using defined equation 
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Abstract of DE1 01 1 2096 

The imaging method has measured data used 
for providing an image, with determination of 
the coordinates of at least one image point at 
which the concentration of at least one matter 
has an extreme value. The time variation in the 
matter concentration for several image 
elements of a first image and a subsequent 
image is determined via a defined equation in 
terms of a Nabla operator and a position- 
dependent diffusion tensor. 
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Prufungsantrag gem. § 44 PatG ist gestellt 

(54) Bildgebungsverfahren und Verfahren zur Rekonstruktion neuronaler Konnektivitat 

Die Erfindung betnfft ein Verfahren, bei dem a us Mess- 
daten wenigstens ein Bild erzeugt wird. 
Erfindungsgemaft wird dieses Verfahren so durchgefuhrt, 
dass Koordinaten wenigstens eines Punktes ermitteit 
werden, fur den die Konzentration von wenigstens einem 
Stoff einen extremal en Wert aufweist, und dass eine zeit- 
liche Veranderung der Konzentration fur mehrere Ele- 
tnente eines ersten BNdes und eines weiteren Bildes er- 
mittelt wird, wofaei eine zeitliche Veranderung der Kon- 
zentration C des Stoffes im Wesentiichen gema£ der For- 
mel 

dC/dt = V > (gVC) 

ermitteft wird r wobei V einen Nablaoperator {in 3-D: 
(d/dx f d/dy, d/dz)^ ) bezeichnet und wobei £ einen 
ortsabhangigen Diffusionstensor bezeichnet. 
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Beschreibung 

[0001] Die Erfmdung betrifft cin Bildgebungs verfahren, bei dem aus Messdaten wenigstens ein Bild erzeugt wird, 
[0002] Bin galtungsgemaBes Verfahren wird in dem Konferenzbeitrag Tracking Fibre Tracts Using Fast Marching, G, 
5 J. M. Parker, Proc. Intl. Soc. Mag. Reson, Med. 8, 85 (2000) beschrieben, Bei dem Verfahren werden durch Losen einer 
Eikonal-Gleichung und durch sukzessive Anderung einer Frondinie wahrscheinliche Verkniipfungen zwischen einem 
Ausgangspunkt und einem oder mehreren Ziclpunkten ermiftelt. 

[0003] Die Wachstumsgeschwindigkcit an einem Punkt der Front ist. proportional zurn Reurag des Skalarprodukts zwi- 
schen dem zu diesem Punkt. gehorenden Front-Normalcn vektor und dem dominanten Eigenvektor des Punktes. Die Zeii, 
to die es dauert, bis die Front einen Punkt Qberqucrt hat (filr eine Anzahl von Bildpunkfen zusammengefasst in einer Time- 
of-Arrival-Map), wird als MaB fur die Konnekti vital zwischen diesem und dem Ursprungs-Punkt. betrachtel. Dieses Ver- 
fahren hat den Vorteil, dass es moglich ist, Verzweigungen zu berucksichtigen. 

[0004] Zur Bestimmung der eigentlichen Trajektorien mittels der Tmie-of- Arrival-Maps erfolgt die Berechnung von 
Gradienten fiir eine Rethe von aufeinanderfolgenden Punkten, die sich jeweils durch Extrapolation entlang des zurn vor- 
ts angegangenem Punkt gehorigen Gradienten ergeben. 

[0005] Dieses und andere Bildgebungsverfahren werten ortsaufgelosle Tnformationen liber die ("^elbst."-)DifT , usion 
wenigstens eines in einer Probe enthaltenen StolTes aus. Ortsaufgeldste Diffusionsdaten werden dabei mil. einem geeig- 
neten Mess verfahren, insbesondere mil: Magnet-Resonanz-Bildgebung, erhaJten. 

[0006] \Selbsr-Diffusion stellt einen Spezialfall der Diffusion dar: Diese, die eine Konsequenz der thermischen Be- 
20 wegung von Atomen, Molekulen u'nd Partikeln (Brown'sche Molekularbewegung) darstellt, fiibrt zu einem Material- 
transport von Gebieten hoher zu Gebieten niedriger Konzentration. Selbst-Diffusion Hegt vor, wenn Konzentration sgra- 
dienten nicht vorhanden bzw. schon ausgeglichen sind. 

[0007] Das Verfahren von Parker und die nachslehend gewurdigten Verfahren zur Bestimmung der (neuronaien) Faser- 
bahnen ( M Fibertracking TI ) verwenden (3 x 3-)Tensoren und vcrschiedene Verfahrensschritte zur Beriicksichtigung von 
25 Diffusions- Anisotropic. 

[0008] Hng verwandt mk dem Verfahren von Sethian sind die Verfahren von Tsitsikiis, dargeiegt in dem Artikel Effi- 
cient Algorithms for Globally Optimal Trajectories: John N, Tsitsikiis, IEEE Trans ac. Aulom. Conlr., 40/9), 1528-38 
(1995). Die dort angegebenen Algorithmen sind nur fiir skalare Fcldcr geeignel, nicht. aber fur Vektor- oder Tensorfclder, 
da fiir diese der Beitrag eines Voxels zurn Pfadintegrai von der lokalen Richtung tier Trajektorie und sotnil. vom zuriick- 

30 iiegenden Trajektorien verlauf abhinge. Das Verfahren, im Folgenden als skalares Tracking- Verfahren bezeichnet, ermit- 
telt minimale Werte fur den Aufwand (Kosten, Energie etc.), der erforderiich ist, urn von einem inneren Punkt an einen 
Randpunkt zu gelangen, Aus diesen Tnformationen konnen dann ciie eigenUichen Trajektorien berechnet werden. 
[0009] Der Artikel Tracking neuronal liber pathways in the living human brain: Thomas H. Conturo, Nicolas F. Lori 
Thomas S. Cult, Erbil Akbudak, Abraham Z. Snyder, Jushua S. Shimony, Robert C. McXinslry, Harold Button and Mar- 

35 cus E. Raichie, Proc, Nail. Acad. Sci. USA, 96, 10422 -10427 (1999) beschreibt ein Verfahren, bet dem fiir jedes Voxel 
inncrhalb einer Startregion mit der Berechnung einer Trajektorie begonnen wird. Entlang der Richtung des "dominanten 
Eigenvektors erfolgt eine lineare Extrapolation uber eine konslanle Lange, die der Halfte der Voxel- Kan ten lange ent« 
spricht [als "dominanler Eigenvektor" eines Tensors wird der zum groBlen Eigenwert gehorige Eigenvektor bezeichnelj. 
Es erfolgt ein Abbruch der Trajeklorienberechnung, wenn die Anisotropie eine bestimmte Grenze unterschreket, das 

40 hetRt , wenn der gewahlte Anisotropie-Index unter einen vorgegebenen Schwellwert fallt. 

[0010] Der Artikel Tracking liber bundles in diffusion tensor images: C Poupon, J.-R Mangin, M. Pachot-Cluart, F. 
Poupon, L Regis, V. Frouin, L Bloch, D, LeRihan, Neurolmage 7 (4), 5701 (1998) beschreibt ein gattungsgemafies Ver- 
fahren, bei dem uberpriift wird, ob zwei benachbarte Voxel zum selben Faserbundel gehoren. Hierzu wird sowohl der 
Winkel zwischen den dominanten Eigenvektoren zweier aufeinanderfolgender Voxel ausgewertel, als audi der Winkei 

45 zwischen zwei aufeinanderfolgenden Voxel verbindungsgeraden. 

[0011] Der Artikel Three- Dimensional Tracking of Axonai Projections in the Brain by Magnetic Resonance Imaging: 
Susumu Mori, Barbara J. Grain, V. P. Chacko and Peter D. M. von Zijl, Ann. Neurol., 45 (2), 265 -270 (1 999) beschreibt 
ein unter der Abkiirzung "FACT" (FACT: Fiber Assignment by Continuous Tracking) bekanntes Verfahren, bei dem eine 
Trajektorie innerhalb eines Voxels durch eine Gerade beschrieben wird, dessen Richtung durch den dominanten Eigen- 

50 vektor gegeben ist, Der Starfpunkt der Geraden des nachslen Voxels definiert sich durch einen Schnitt der Geraden des 
aktueilen Voxels mit der Oberflache dieses Voxels. Das Fiber-Tracking fur einen Fascrstrang wird abgebrochen, wenn 
der Miltelwert der Skalarprodukte zwischen dem dominanten Eigenvektor des aktueilen Voxels mit den normierten, do- 
minanten Eigenvektoren der drei nachsten Voxel unter einen gegebenen Schwellwert fallt (zum Beispiel 0.75). 
[0012] Eine Anwendung des obigen Verfahrens zur Bestimmung des Verlaufs von Faserbahnen im Raltenhirn wird 

55 z, B. beschrieben in "In Vivo Three-Dimensional Reconstruction of Rat Brain Axonai Projections by Diffusion Tensor 
Imaging: Rong Xue, Peter C. M. von Zijl, Barbara J. Grain, Meiyappan Solaiyappan and Susumu Mori, Magn, Reson. 
Med. 42 (6), 1123-27 (1999) 1 '. 

[0013] In dem Artikel. Non-Invasive Assessment of Axonai Fiber Connectivity in the Human Brain Via Diffusion Ten- 
sor MRI: K. Derek Jones, Andrew Simmons, Steve C. R, William und Mark A, Horsfield, Magn. Reson. Med., 42 (1), 
60 37-41 (1999) ist ein Segmentierungs verfahren dargestellt, das Cluster von miteinander verbundenen Voxeln ermittclt. 
Bet dieser Mcthode wird angenommen, dass der Faservcrlauf lokal zwischen zwei Voxein einem Kreisbogen folgt. 
[0014] Fur die Visualisierung von aero- oder hydrodynamischen Slromungsfeldern findet folgende gewohnliche Diffe- 
rentia l-GIeichung Anwendung (Streamline-Ansatz): 

65 dx(x 9 T)fdi = v('x, r). 

[0015] Eine Verwendung dieser Gleichung zur Visualisierung des Blutflusses ist. aus dem Artikel Visualizing Blood 
Flow Patterns Using Streamlines, Arrows, and Particle Paths: Michael II. Buonocore, Magn. Reson. Med., 40, 210 226 
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(1998) bekannt. 

[0016] Eine auf obiger Gleichung basierendes DT4'ibertrac king- Verfahren ist in dem Artikel In Vivo Fiber Tractogra- 
phy Using DT-MRI Data: Peter J. Basser, Sinisa Pajevic, Carlo Pierpaoli, Jeffrey Duda und Akram ALdroubi, Mag. Re- 
son, Med., 44, 625 -632 (2000) beschrieben worden. Hierbei wird als Geschwindigkeii. v der dominante Eigen vektor des 
Diffusion slensors gewahlt. 5 
[0017] Bin Merkmal dieses Verfahrens ist eine hone Empfindiichkeit gegeniiber lokalen, zum Beispiel durch Rauschen 
bedingten Storungcn der Richtung des dominanten Eigenvcktors. Hieraus ergibt sich der in der genannten Arbeit be- 
schriebcne Nachteil, dass kleine Variationen der Startposition zu erhebiichen Unterschieden in den berechneten Trajek- 
torien fuhren konnen, 

[001H] Der Grand fur diese Instability, die auch anderen bekanntcn Tracking- Verfahren zu eigen ist, besteht darin, 10 
dass die Berechnungen fur die einzelnen Teile einer Trajektorie miteinander verkettef sind und fur die Berechnung eines 
Teilstucks nur die Information eines Voxels oder einer sehr kleinen Umgebung (zum Beispiel einer 1-VoxeFUmgebung) 
verwendet wird, 

[0019] Das Verfahren von Basser et. at, und ahniiche Verfahren setzen die Stetigkeit von benachbarlen Eigenvektorcn 
voraus, die jedoch aus mindestens zwei Grunden in vielen Fallen nicht gegeben ist: 1. Eigenvektorcn des Tensors sind is 
axiale und keine polare Vektoren, d. h. die Vektoren besttzen keine ausgezeichnete Richtung, wodurch Vorzeichenwech- 
sel zwischen benachbarten Eigenvektorcn eines Paserbundels auflreten konnen. 2. Wenn sich aufgrund von Rauschen die 
Reihenfolge der beiden grMten Eigenwerte andert, hat dies eine Vertauschung der beiden zugehorigen Eigen vektoren 
zur Folge, Da die Eigenvektorcn zueinander orthogonal sind, konnen somit auch Richtungsanderungen urn 90° zwischen 
zwei benachbarten Eigenvektorcn einer Trajektorie auftreten, insbesondere bei hinreichend kleiner Anisotropic. 20 
[0020] Auf beide Aspekte, und insbesondere den zweiten, wurde eingegangen in dem Artikel Statistical Artifacts in 
Diffusion Tensor MRI (DT-MRI) Caused by Background Noise: Peter L Basser und Sinisa Pajevic, Magn. Reson. Med, 
44,41 50(2000), 

[0021] Aus dern Artikel Regularization of Diffusion-Based Direction Maps for the Tracking of Brain White Matter 
Fascicles: C. Poupon, C A. Clark, V. Frouin, J. Regis, I. Bloch, D. Ix Biban unci J.-F. Mangin, Neurol mage 12, 1 84-1 95 25 
(2000) ist.ein Verfahren bekannt, bei dem eine Berechnung von neuen Vektorfeldern durch lokale Miltelung von Eigen- 
vektorcn erfolgt. Wie sich jedoch zeigt, hangen die ftir diese Vektorfelder berechneten Trajektorien erheblich vom Grad 
der Mittelung ab. 

[0022] Bin Nachteil der oben genannten Verfahren ist, dass die durch den Diffusions-Tensor zur Veriugung gcstcllte In™ 
formation nur zum Teit genutzt wird. Bin Ansatz, der fur die Trajektorien-Berechnung die gesamte Information des Dif- 30 
fusion st en sors nutzt, ist in dem Konferenzbeitrag Axon Tractography with Tensorlines: Mariana Lazar, David Weinstein, 
Khadcr Hasan, Andrew L. Alexander, Proc. Intl. Sac. Mag. Reson. Med., 8, 482 (2000) dargestelk. Zunachst wird mit 
dem Diffusionstensor eine neue Richtung defmietl, 

Yout = C ' &n , (1) 35 

wobei Vj tl die ursprungliche Richtung der Trajektorie angibt, Aus v OER wird dann wie folgt ein neuer Richtungsvekior be- 
rechnet: 

}£pro P = d-a)- & n + a ■ ^ (2) 4 o 



[0023] Hierbei stcltt a einen Gewich lungs faktor dat; fur den ein Wertebereich von 0.1 bis 0.6 vorgeschlagen wurde. 
Dieser Vektor ersetzt den dominanten Eigen vektor beim konven Lionel len Streamline- Ansatz (siehe oben). 
[0024] Ein weiteres Verfahren wird in dem Artikel A Path Integral Approach to While Matter Tractography; David vS. 
Tuch, John W, Belli veau und Van J. Wedeen, Proc. Intl. Soc. Magn. Reson. Med., 8, 791 (2000) beschrieben. Es basiert 45 
auf einem Pf ad-Integral- Ansatz, bei dem ein Trajektorien- Verlauf zwischen zwei Punkten simultan optimiert wird. Fiir 
die Anwendung des Verfahrens ist die Bestimmung mehrerer PotentiaL-Funktionen erforderlich, Aufgrund der damit ver- 
bundenen Schwierigkeiten sind aber offenbar bis jetzt noch keine Impiementierungen dieses Fiber-Tracking- Ansatzes 
realisiert worden. 

[0025] Der Erfindung liegt. die Aufgabe zugrunde, ein gattungsgemaBes Verfahren bereitzustellen, durch das das Vor- 50 
handensein faserartiger Strukturen in einer Probe nachgewiesen werden kann, mit dem sich feststellen lasst, ob zwei be- 
liebige Punkte innerhalb der Probe durch die Faserstruktur miteinander verbunden sind und mit dem der Verlauf von fa- 
serartigen Strukturen zuveriassig erfasst und dargestelk werden kanm 

[0026] ErfindungsgemaG wird diese Aufgabe dadurch gelost, dass ein gattungsgemaBes Bildgebungsverfahren gemaB 
den Merkmalen des Paten tans pruchs 1 durchgefuhrt. wird, 55 
[0027] Die Erfindung sieht insbesondere vor, das Bildgebungsverfahren so durchzu ftihren, dass Diffusion sdaten auf- 
genommen werden, dass die Koordinaten wenigstens eines Punktes in einer Probe bestimmt werden, der zur Festlegung 
einer Anfangs-Konzentrationsverteilung wenigstens eines Diffusionsmediums mit wenigstens einem extremalen Wert 
verwendet wird t dass cine zeitliche Veranderung der Konzentration fur mebrere Blemente eines ersten Biides ermittelt 
und als ein weiteres Bild erzeugt wird, wobei eine zeitliche Veranderung der Konzentration C des Diffusionsmediums im 60 
Wesentlichen gemaB der Formel 

dc/dt = v . (gvo 

bestimmt wird. Hierbei bezeichnet V den Nab la- Operator (dfdx, d/dy, dld't)' 1 und D einen (ortsabhangigen) Diffusions- 65 
tensor. 

[0028] Eine bevorzugte Ausfuhrungsform der Erfindung sieht vor, dass Linien oder Fiachen konstanter Konzentration 
beschrieben werden und dass in jedem Zeitsc brill mindestens ein Punkt. auf diesen Einien oder Fiachen zur Bestimmung 
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einer Trajektorie ausgcwahlt wird. 

[0029] Bs ist. besonders zweckmaftig, das Verfahren so durchzufiihren, dass die Gleichung 
dC/dt = V • mVC) 

5 

fiir verschiedene Zeiten gelost. wird, wobei fur wenigstens eine crste Zeil \q eine Ausgangskonzentration (Co) zugrunde 
gelegt wird, die ein ausgepragtes Maximum aufweist. 

[0030] Die Verteilungsfunktion ist so zu wahlen, dass sie in der Umgebung des Maximums zumindest. naherungsweise 
eine Radial-Symmetric aufweist, Ais besonders zweckmaBig erweist es sich, wenn die Verteilungsfunktion im Wesentli- 
10 chen der Formel 

Co - exp(- ((x - x 0 ) 2 + (y - y 0 ) 2 + (z - z<,) 2 )/(2<J 2 )) 

entspricht, in der C die Breite der Verteilungsfunktion angibt, die zweckmaBigerweise kieiner oder gleich der Langendi- 
15 mension eines Voxels gewahlt wird und bei der x 0) y 0 und z 0 die kartesischen Koordinaten des (Peak-)Maximums ange- 
ben, das an den Punkt, fur den die Konnektivital untersucht werden soil, oder an einen benachbarten Punkt gelegt wird. 
[0031] Vorzugsweise wird das Verfahren so durchgefuhrt, dass die Komponenten der Diffusionstensoren, die sich in 
Bereichen mil nahezu freier Diffusion, also groBen Werten der Spur des Tensors befinden, auf einen sehr kleinen Wert 
(z. B. null) gesetzt werden, 

20 [0032] Fernerrist es zweckmaBig, dass die Komponenten der Diffusionstensoren, die sich in Bereichen niedriger Dif- 
fusion-Anisotropie befinden, auf einen sehr kleinen Wert (z. B. null) gesetzt gesetzt werden. Vorzugsweise fmdet der 
RA-Index (RA: Relative Anisolropy) Anwendung, der definiert ist als RA = 1//3 * ((X { - k ba( .) 2 + (X 2 X hdX ) 2 + (X 3 ~ 
X h&l ) 2 ) U2 IX b ^ mil hx- = (h + X 2 + X 3 )B. 

[0033] Auch ist es vorteilhaft, dass die Komponenten der Diffusionstensoren mil einem (skalaren) Diffusions- A niso- 
25 tropic-Index gewichtet werden. 

[0034] Femer ist es vorteilhaft, dass Linien oder Flachen konst.an1.er Konzefflralion (Isolinien, -flachen) besdmmt und 
Abweichungen von einer Kreis- bzw, Kugelform ermiuelt werden, wobei sich fur die nachfolgenden Verfahrensschrilte 
eine Glattung der Tsolinien bzw, * flachen als zweckmaBig erweisen kann, 

[0035] AuGerdem ist es vorteilhaft, dass in jedem Zcilschritt die Konzentrationswerte der Isolinien bzw. -flachen so ge- 
30 wahlt werden, dass sie einem konstanten Bruchleil (z, B. 0,607, siche unten) des momentanen Maximalwertes der Kon- 
zentradonsfunktion oder des aktueilen Wertes der Konzenlralionsfunktion am Ort des anfangltchen Konzentrationsma- 
ximums entsprechea 

[0036] Die Messung des DifTusionstensors kann auf verschiedene Weisc erfolgen. Fiir die Untcrsuchung von Diffusi- 
onsprozessen in biologischen Systenten, sind MR-Diffusions-Messungen am bcsten gceignet, da sie zur Zeit die einzige 
35 Moglichkeit darstellen, nichl-invasiv die Bewegung von Molekulen Liber Distanzen zu verfolgen, die (in einem slalisti- 
schen Sinn) mit Zellabmessungen vergleichbar sind. 

[0037] Fur Messungen der Diffusion mittels Magnetresonanz (MR) wird dem Hauptmagnetfeld ein schwacheres Ma- 
gnelfcld uberlagert, dessen Starke im Wescnilichen linear Liber den Oil variiert (Diffusionsgradient). Die Bewegung von 
Kernspins in diesem Gradienlenfeld fiihrt zu einer Abnahme des gemessenen MR-Signals, aus der die Starke der Diffu- 
40 sion berechnet wird, 

[0038] Die Messung der Diffusionsprozesse mittels kemmagnetischer Resonanz erfolgt. in einer besonders geeigneten 
Weise durch BUdgebimgssequenzen, die auf dem im Artikel Spin Echoes: E. L. TIahn, Physical Review, 80/4) (1950) 
oder dem in Spin Diffusion Measurements: Spin Echoes in the Presence of a Time-Dependent Field Gradient: E. O. 
Stejskal und J. E. Tanner, J. Chem. Phys. 42 (1), 288- 292 (1965) beschriebenen Verfahren basieren. Die Messung der 
45 Diffusionstensoren erfolgt bevorzugt mittels eines Bildgebungsverfahrens, das sich an dem in dem Arlikel Use of Spin 
Echoes in a Pulsed Magnetic-Field Gradient to Study Anisotropic, Restricted Diffusion and Flow, E. 0, Stejskal, J. 
Chem. Phys, 43 (10), 3597-3603 (1965) beschriebenen Verfahren orientiert. In dieser Arbeit werden Ausdriicke fur die 
Abhangigkeitder Signal-Amplitude von der gewahlten Gradientenrichtung angegeben. 

[0039] Mit. MR-Diffusionsmessungen am lebcnden Gewebe wird ublicherweise die Seibstdiffusion von Wasser unter- 
50 sucht, die Ruckschlusse auf verschiedene physiologisch relevant© Gewehe-Eigenschaften und -Parameter eriaubt, wie 
zum Beispiel Temperatur, Permeabilif at der Zellmembranen, AusmaB der Diffusionshinderung und Anisotropic. Letztere 
tritt insbesondere in und um Nervenfasem auf: Die Diffusion entlang der Nervenfasern erfolgt deutlich schneiler als 
senkrecht zu deren Richtung. 

[0040] MR-DilTusionsmessungen besitzen verschiedene klinische Anwendungen, Als Einsatzgebiete sind neben der 
55 Schlaganfalldiagnostik die Untcrsuchung von Tumoren und Multipler-Sklerose zu nennem 

[0041] Fiir die Bestimmung der Komponenten des DifTusionstensors werden aufeinanderfolgende Di (fusion s-Mcssun- 
gen durchgefuhrt, bei denen von einer zur nachsten die Richlung des DifTusionsgradienten variiert wird. 
[0042] Das erfindungsgemafle Verfahren wird vorzugsweise fur die Rekonslruktion neuronaler Konnektivital aus DT- 
MRI-Daten eingesetzt, wobei mit neuronaler Konnektivitat die Vernelzung verschiedener Cjehirnareale durch Nerven- 
60 bahnen (in der weiBen Hirnsubstanz) gemeint ist. 

[0043] Weitere mogliche Anwendungen des Verfahrens bestehen in der Untcrsuchung von anderen nat urlichen Fasem, 
wie zum Beispiel Muskelfasern, ktinstlich erzeugten Fasern und hieraus gewonnenem Gewebe sowie von kiinstlichen 
oder natiirlichen porosen Medien (worunter zum Beispiel auch Gefafte und Kapillare in derischen und pflanzlichen Or- 
ganismen verstanden werden sollen), 
65 [0044] Die Erfindung beinhaltet den Einsatz einer Simulation einer Ausbreitung von virtuellen Konzentrationsvertei- 
iungen zur btldgebenden Darstellung von Konzentrationen und/oder Konzcntrationsanderungen oder zur bildgebenden 
Darstellung von daraus abgeleiteten GroBen. 

[6045] Es ist zweckmaBig, Konzcntrationsanderungen ausgehend von einer Startregion durch ein Wachstum von ge- 
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scblossenen Frontlinien zu erfassen unci darzustellen und aus einer Reihe von aufeinanderfolgenden Frontlinien (bzw. - 
flachen) Punkte auszuwahlen und zwischen diesen zu interpolieren, 

[0046] Weitere Vorteile, Besonderheiten und zweckmaBige Weiterbildungen der Erfindung ergeben sich aus den Un- 

tcranspruchen und der nachfolgenden Darstellung bevorzugter Ausfuhrungsbeispiele anhand der Zeichnungen. 

[0047] Von den Zeichnungen zeigt 5 

[0048] Fig, 1 eineDiagramm-Darstellung der Stejskal-Tanner-Sequenz, 

[0049] Fig. 2 ein EPI-Sequenz-Diagramm mil zugehorigcm k-Raum-Abtast-Scherna, 

[00S0] Fig. 3 die PuLssequenz zu einer Singlc-Shot-STEAM-Sequenz, 

[0051] Fig. 4 im obcren Teilbild eine Darstellung eines RA-Index-Map, das fur jedes Voxel der aufgenomrnencn Ge- 
hirnschicht den zugehorigen RA-Index afs Xntensitatswert darsteilt (Voxel mit RA < 0.1 sind schwarz dargestelit) to 
und irn unteren Teilbild eine Moment aufnahme der Diffusion sausbreitung eines anfanglichen Konzentration speaks. 
[0052] Fig* 5 Eine Ausbreitung eines Stoffes in einer Faser, 

[0053] Fig. 6 eine Aufspallung eines Faserbundels und Isolinien von Konzentrationsverteilungen innerhalb der Faser 
sowie 

[0054] Fig* 7 in funf Teiibildern A, B, C, D, E Isoiinien bei Diffusions vorgangen analog zu den Fig* 6 dargestellten unci is 
zugehorige Kriimmungswerlc. 

[0055] Besonders zweckmaBige Ausfiihrungsfonnen der Erfindung sehen fur die Durchfuhrung der Diffusionsmessun- 
gen die Verwendung von Single-Shot- BiLdgebungssequenzen von Es ist insbesondere zweckma'Big, eine Sequenz zu ver- 
wenden, die aus dem in Fig. 1 dargestellten DifTusionsgewichtungs-( M Slejskal-Tanner ,, -)rVlodul und dem nachgeschalte- 
ten Bildgebungsmodui aus Fig. 2 besteht oder eine Fig. 3 entsprechende Sequenz zu verwenden (Single-Shot-STEAM). 20 
Zusalzlich zu den beiden Gradienlen-"Pulsen" ist ein statiscber Hintergrund-Gradienl: dargestelit, der z. B. durch Suszep- 
tibiiitatssprunge zwischen zwei verschiedenen Geweben verursacht sein kann. 

[0056] Bei Kombination des Stejskal-Tanner-Moduls mil der in Fig. 2 (iinks) dargestellten Gradienten-Echo-Planar- 
Bildgebungssequenz, wie sie z, B. fur funktionale MR-Gehirn-Bildgebung verwendet wird, entfallt der eingezeichnete 
HF-Anregungspuls. Das zur Sequenz gehorige k-Raum-Abtasl .schema ist in Fig, 2 (rectus) dargestelit, 25 
[0057] Die STEAM-Sequenz (Fig* 3) besteht. aus einem HF-Anregungspuls, dem DitTusionsgewichtungs-Modul, ei- 
neni 90°-Puls, mit dem die transversale wieder in eine longiludinale Magnetisierung verwandelt wird, und einer Serie 
von IIF-Anregungspulsen mit kleinem Flip-Winkel. Diese Sequenz, die gegeniiber einer EPI-Single-Shot-Sequenz ein 
um ca. 50% reduziertes SNR aufweist, findet: vorzugsweise bei BO-Feldstarken von mehr als 1.5 T Verwendung, Der 
(jrund dafiir ist, dass sich mit. zunehmender BO-Feldstarke die transversale Reiaxalionszeit T2* verkurzt, wodurch EPf 30 
aufgrund der steigenden Anforderungen an die Gradienten performance zunehmend unpraktikabel wird. 
[0058] Eine zweckmaBige Ausfuhrung der Erfindung sieht vor, die Ausbreitung eines virtue] len Konzentration speaks 
ftir einen gemessenen Datensatz von (ortsabhangigen) DitTusionskoeffizienten oder Diffusion stensoren zu berechnen. 
Dieses erfoigt durch das Iiteen der part iel len Different! al~Gieichung 



dC/dt = DV 2 C, (3) 
irn Fall von isotroper Diffusion, oder fur anisotrope Diffusion durch das Losen von 

dC/dt = V • (gvc) . (4) 

[0059] In Komponenten-Schreibweise lautel diese Gleichung: 

dC/dt - Z 3 iij . a B 1 j ^C/ax^Xj (5) 



Go - exp(-«x ~ x 0 ) 2 + (y - y 0 ) 2 + (z - ^) 2 )/(2g 2 )) 



35 



40 



45 



[0060] Jede der drei Gleichungcn (3-5) wird nachfolgend a Is T, Diffusionsgleichung" bezeichnet. 

[0061] Zur vollstandigen Definition von DilTerenliaigleichungen gehort zusat/Jich die Angabe der Anfangs- und Rand- 

bedingungem Erst ere werden fur das gesamte Berechnungsvolumen angegeben, letztere nur fur dessen Oberflache, 

[0062] Die Tensordaten D konnen sich sowohl auf ein TeilvoLumen des untersuchten Objekis a!s auch das untersuchte 

Objckt einschlieBlich des ihn umgebenden Raums (in dem sich die DifFusionskoeffizienten als Null annehmen iassen) 50 

bezichen. 

[0063] Iris erweist sich als zweckmaBig, Neumann- Randbedingungen und speziell den Fluss durch die Oberflache als 
null anzunehmen. Wenn jedoch die Diffusionskoeffizienten/-tensoren im Randbereich null sind, ist es unerheblich, ob 
eine Neumann- oder eine Dirichiet-Randbedingung gewahlt wird. 

[0064] Die Anfangsbedingung ist gegeben durch eine Ausgangskonzentration (zum Zeilpunkt to), die ein ausgepragtes 55 
Maximum aufweist. Die Verteilungsfunktion sei in der Uingebung des Maximums zumindest naherungsweise radial- 
syntmetrisch. Als besonders zweckmaBig erweist es sich, wenn die Verteilungsfunktion im Wesentlichen der Formel 



60 



entspricht, in der a die Breite der Verteilungsfunktion angibt. und x 0 , yo sowie z 0 die kartesischen Koordinaten des 
(Peak-)Maximums bezeichnen. Die Breite a wird zweckmaRigerweise kleiner oder gleich der Langendimension eines 
Voxels gewahlt. Das Peak-Maximum wird an den Punkt gelegt, dessen Konnektivitat untersucht werden soil (oder in des- 
sen Nahe), 

[0065] Losungen der DifTusionsgleichungen iassen sich, wie allgemein bei lincaren partiellen DilTerentialgleichungen, 65 
durch eine Diskretisierung und durch I^xsung der sich daraus ergebenden linearen Crleichungssysteme erhallen. Hierbei 
konnen Finite-Differenzen-Methoden oder Finite-Elemente-Methoden (FEM) verwendet werden. Die Generierung der 
(jitter, auf denen die Berechnungen durchgefiihrt werden, kann zum Beispiel mitteles eines Deiaunay-Algorilhmus er- 
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folgen. 

[0066] Zur Minimierung der Rechenzeit ist es zweckmaBig, das Gitter so zu wahlen, dass es in Bereichen mil einer 
starken Variabilitat des Diffusionskoeffizienten eine groBere Knotendichte aufweist als in Gebieten mil einem relaiiv ho- 
rn ogenen Diffusionskoeffizienten (/-tensor). 
5 [0067] Zur Losung der DifYusionsgleichung konnen verschiedene Programme verwendet werden, wobei der Fachitiann 
vorzugsweise bekannle Program mpakete einsetzt, wie beispielsweise "difTpack 11 (http://www.numerical-objecls.com/) 
oder "femlab" (hup:/ / www.femlab.com). Bine umfassendere AufUstung von Programmpaketen fur das Losen pariieller 
Differentialgleichungen findet sich unter: Internet Finite Element Resources: Young, Roger and MacPhcdran, Ian, 
<http:/ / www,engr.usaskxa/-inacphed/finite/fe resources/fe_rcsources.htini> (2001). 

io [0068] Urn zu vermeiden, dass beim Fiber-Tracking jenes Gewebe rait cinbezogen wird, das die Faserbtindel umgibt, 
ist es zweckmaBig, die Diffusions tensoren fur solche Gebiete auf null zu setzen, in denen die Diffusionsanisolropie (zum 
Beispiei beschrieben durch den Relative- A nisotropy-Index) einen gcgebenen ScbweUwert unterscbreitet. Alternativ 
konnen diese Gebiete durch das Einfiihren innerer Grenzflachen ausgesparl werden. Der "Relative Anisotropy Index", 
RA, ist dcfinicit als RA = 1//3 - ((k { - X ba! .) 2 + (h - X h ^ 2 + (X 3 -A^) 2 ) 172 /^ wobei X bar - (K { + X 2 + X 3 )/3 unci X, die 

1 5 drei Bigenwerte des DifFusionstensors bezeichnet. 

[0069] Anstelle einen Schwellwerl zu setzen, kann es vorteilhaft sein, die Komponenlen der Diffusionstensoren rait ei- 
nem (skalaren) Diffusions- A nisofropie-Index zu gewichten, wofiir vorzugsweise ebenfalls der RA-Index verwendet wirci 
und wobei eine vorherige Norraierung der Diffusions- Tensoren zweckmaBig ist, die vorzugsweise dadurch erfolgl, dass 
dieBlemente eines Tensors durch dessen Spur dividiert werden. 

20 [0070] EiB-Beispie! fur eine zweidimerisionale Ausbreitungsberechnung fur einen plariaren, das heiBt zweidimensio- 
naien Diffusionsdaiensatz von einem menschlichen Gehirn ist in Fig. 4 dargestelit. In der oberen Bildhalfte ist das RA- 
Index-Map gezeigt (bei dent zur Hervorhebung der weiBen Hirnsubstanz ein Schwellwert gesetzt wurde) in der unteren 
eine Mornentaufnahme der Diffusion sausbreitung des anfanglichen Konzent rations-Peaks, dessen Position durch den 
schwarzen Punkt in der oberen Bildhalfte dargesteih ist (ira Corpus Callosiim). 

25 [0071] Bine weitere Moglichkeit fur das Losen der DitTerenlial-Gleichung besteht in einer Nachbildung des Simulati- 
onsproblems durch geeignete physikalische Model [-Systeme, zum Beispiei durch analog-elektronische Schaltkreise, 
[0072] Fur das eigentiiche Fiber-Tracking wird dann die Ausbreit.ung von Konzentrations-Isolinien oder -Isoflachen 
wahrend des siraulierten Diffusionsprozesses verfolgt (siehe Fig. 5). Zur Ermittlung der Lsoiinien, beziehungsweisc Iso- 
flachen, ist in jedem Simulationsschritt ein neuer, gegeniiber dem vorangegangenen Schritt kleinerer Konzent ratio nswert 

30 festzuiegen (Konzentrations-Isolinien: Verbindung von Punkten mil gleicher Konzentration), Dieser Wert kann bei- 
spielsweise so gewahlt werden, dass er iramer einen konstanten Bruchteil der augenblicklichen Maximal- Konzent rat ion 
(bzw. der Konzentration an! Ort des Anfangspeaks) betragt. 

[0073] Bine weitere Moglichkeit zur Festlegung der Isolinien/flachen besteht darin, die Punkte zu bestiramen, bei der 
die Krummung der Konzentrationsfunktion verschwindet oder minimal wird (Kriimmungs-Isoiinien). In einer Dimen- 
35 sion entspricht die Krummung der zweiten Ableitung der Konzenlraiionsfunktion bezuglich des Orts. In zwei und mehr 
Dimcnsionen wird die Krummung veraiigemeincrl zu k c = (grad(C)) T H(C)(grad(C)), wobei II die Hesse-Matrix zur 
Konzcnlrationsfunktion C bezeichnet (Hy = £ 2 Xjdxj). Im zweidiraensionalen Fall gill also: 

a fl r = (dC/dx i , dadx 2 )K(dOdx i ) 2 + (dcm 2 ) 2 y 1/2 

40 

h r= {(eKifdx^ d 2 adx\dx 2 \ (e 2 adx x dx % a 2 c/£x 2 2 )) 

k c ~ aV R ao, 

45 wobei II in Zeil en notation angegeben ist. 

[0074] Im Fall einer normierten Gaussfunktion (p) in einer Dimension und bei isotroper Diffusion ergeben sich aus der 
Bedingung d 2 p/dx 2 = 0 die x-Koordinaten x = ±/2j/t, Hieraus folgt, dass das Kruramungskriterium und das Ampiit.uden- 
kriterium ubereinstimmen, wenn fur den Amplitudcnschweliweri der Wert. exp( - 1/2) w 0.6065 gewahlt wird. 
[O075] Als MaB fur die Konncktivhat des Ausgangspunktes zu einem anderen Punkt kann wie bei der Fast-Marching- 

50 Methode (10) die Zeit definiert werden, die die Isolinie, beziehungsweisc Isoflache, bis zum Hrreichen dieses Punktes be- 
notigt (t adv?l ). 

[0076] Diese Daten werden in einem Time-of- Arrival-Map (bzw. -Array) zusammengefasst. Ausgehend von einem 
Startptmkt werden fur eine Reihe von aufeinanderfolgenden Punkten Gradienten der Time-of- Arrival-Maps bereebnet, 
wobei der jeweils nachste Punkt als ein Punkt der Geraden gewahlt wird, dessen Aufpunkt dem aktuellen Punkt ent- 
55 spricht und dessen Richtungsvektor durch den CJradienten fur den aktuellen Punkt gegeben ist, Auf diese Weise ergibt 
sich eine Trajektorie, die den gewahlten Startpunkt mil. dem Ort des Anfangspeaks der Ausbreitungsberechnung verbin- 
det ("back-propagation'^-analog zur Fast-Marc hing-Methode. 

[0077] Bei isotroper Diffusion, d. h. bei gaussformiger Diffusionsausbreitung ist Proportional zum Quadrat des 
Abstandes vom Ort des ursprunglichen Konzentrationsmaximums (Zentrum). Dies fuhrtdazu, dass die Zeit fiir das tjber- 

60 streichen eines Trajektorienstiicks rait zunebmendem Absland vom Zentrum anwachst. Bei der Verwendung der Timc- 
of- Arrival-Maps kann es daher sinnvoll sein, die Zeit Wval durch die effekiive Zeit tarrival r anival - i/t^vai zu ersetzen. 
[0078] Bine weitere Moglichkeit ftir das Aufstellen von Trajektorien besteht darin, von den Lsoiinien/- flachen fur auf- 
einanderfoigende Zeitpunkle jeweils den Punkt auszuwahlen, der den maximalen Absland zum Ursprung aufweist, Zwi- 
schen dtesen Punkten wird dann durch geeignete Funktioncn, wie z. B. (jeraden, interpolicrt. 

65 [0079] Zur Verringerung des Rechenaufwands bei der Bestimmung der Ii3sungen der partieUen Differentialgleichung 
kann es sinnvoll sein, die Ausbreitungsberechnung nichl. fiir den Diffusion stensor-Datensatz als ganzes durchzuftihren, 
sondern fur eine Reihe von uberlappenden Teilvoiumen. Derartige Teilvolumen sind in Fig. 6 beispielhaft dargestelit, 
[0080] Die Definition der Start- und Hndpunkte von Trajektorien konnen zum Beispiei auf der Basis von anatomi- 
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schem Wtssen manuell gewahlt werden odcr an hand von fMRI-Akdvierungskartem 

[0081] Zur Hrkennung etner Aufspaltung der in Fig. 6 dargestellten Faserbiindel erfolgt eine Bewertung dcr Krum- 
mungen fur die Isolinien oder Isoflachen, wobei die Krummungen beispielhaft in Fig* 7 dargestellt sind und auf einc dent 
Fachmann bekannlen Weise berechnel werden [siehe z. B. Taschenbuch der Mathematik, Bronstein, Semendjajew, Mu- 
siol unci Muhlig, Verlag Hard Deutsch, 1993]. 5 
[0082] Pur kartesische Koordinaten wird die Krummung von Isolinien, die durch Funktionen y(x) beschrieben werden, 
im Wesentlichen gemaS der Formel 

k d 2 .v/dx 3 

[l + (dy/dx) 2 f 2 10 

berechnct 

[0083] In Fig. 7 ist fur zwei verschiedene Isolinien (A und D) die Krummung als Funktion von x dargestellt. Fig* 7B 
zeigt die Krummung fur die Halbellipse aus Fig. 7A. Aus der Halbellipse wurde durch Mullipiikarion mil der Funktion 
aus Fig* 7C die in Fig* 7D gezeigle Kurve generietl, die eine idealisierte IsoJinie flir eine DifTusionsausbreiiung an einetn 15 
Verzweigungspunkt darstellL Die zugehorige Krummung sfunkdon ist aus Fig, 7B ersichtlich. 

[0084] Im Fall von Isoflachen, die <lurch Funktionen z(x, y) beschrieben werden, erfolgt fur die Punkte (x, y) die Be- 
rechnung zweier Krummungsradien, wozu eine quadratische Gleichung gelost wird, die im Wesentlichen die foigende 
Form hat: 



(rt - s z )R 2 + h[2pqs (1 -J- p 2 )t • ( 1 4- q>]R + h 4 - 0, 



h-yi + p2 + q2 

[0085] Aus diesen Kriimrmingsradien erfolgt die Bcrechnung einer mittleren Krummung, im Wesentlichen gemaB der 
Formel 

K = V4(l/Ri + MR2). 



20 



mil 

p = d'/Jdx, 

q = d'/Jdy, 25 

r = dhJdx", 

s = dhJd\dy, 

t = BhJBy 2 

und 



30 



35 



[0086] Fm Ailgemeinen konnen die auftretenden Ableitungen nicht anaiytisch sondern nur numerisch berechnet wer- 
den. Verfahren hierfiir sind bekannt aus dem Buch Numerical Recipes in C, 2 ml edition, W t IL Press, S. A. Teukolsky, W. 
T Vettcrling, und B. P. Flannery, Cambridge University Press. 40 
[0087] lis erweist sich als zweckmaBig, zur Fesliegung von Anfangspunkten (mil den Koordinaten xo, yo und zo) fur 
die Diffusions-Simulation skalare Tracking -Verfahren zu vcrwenden, also Verfahren die Trajektorien an hand von skala- 
ren Dif fusions- Anisotropie-Daten bestimmen. Diese Anisotropie-Daten werden aus den DilTusionstensor-Daten z, B, 
durch die Bcrechnung des RA-Index erhaltem 

[0088] ITierbei ist es besonders zweckmaBig, dass das skalare Tracking verfahren so konzipiert ist, dass es das Pfadin- 45 
tegral fiir eine Trajektorie zwischen zwei Punkten minimiert, wobei uber das Inverse eines AnisotropieTndexes (vor- 
zugsweise den RA-Index) integriert wird. Die Optimierung kann mittels der "Fast Marching M -Memode (von Serbian) er- 
folgen, oder mit den daunt niathcmatisch eng verwandten Methoden, die in dem Artikel Efficient: Algorithms for Glo- 
bally Optimal Trajectories: John N. Tsitsiklis, IBBETransac. Autom. Contr,, 40/9), 1528-38 (1995) beschrieben sind. 
[0089] Bin Hauptvorteil des erfmdungsgemaBen Verfahrens besteht darin, dass es auf die Verwendung von Bigenvck- 50 
toren verzichtet. Bei der Verwendung von Figenvektoren gibt es das Problem, dass die Richtung der Figenvektoren nicht 
eindeutig definiert ist und dass die Sortierung der .Figenvektoren entsprechend der GroBe des zugehorigen Eigenwertes 
anfallig ist fiir Rauschen. Abhangig von der numerischen Stabilitat der jeweiiigen Implememierung der Di [fusion s-Si- 
mulation wird die Empfmdlichkeit gegenuber Rauschen auch dadurch reduziert, dass das virtuelle Diffusionsmedium um 
Storungen "herum" diffundiercn kann, Auch zeichnet sich das erfmdungsgemaBe Verfahren dadurch aus, dass die Infor- 55 
mation der Diffusions-Tensoren praktisch vollstandig gcnutzt wird (bis auf Tensorcn* deren Element© auf null gesetzt 
wurden). 

[0090] Bin besonderer Vorteil der Erfindung ist, dass die simulierte Diffusion auch experimentell nicht nachwetsbare 
Diffusionsprozesse aufzeigt. Die mit MR-Messungen zugangiiche Diffusion erstreckt sich maximal uber ca. \ mm, wo- 
hingegen die simulierte Diffusion Langen von mehreren cm uberbriicken kann, 60 
[0091] Das erfmdungsgemaBe Verfahren wird vorzugsweise fur die Rekonstruktion neuronaler Konnekti vital, aus DT- 
MRI-Daten eingesetzt, wobei mit neuronaler Konnektivitat die Vemetzung vcrschiedener Gehirnareale durch Nerven- 
bahnen (in der weiBen Hirnsubstanz) gemeint ist. 

[0092] Weitere moglichc Anwendungen des Verfahrens bestehen in der Untcrsuchung von anderen natiiriichen Fasern, 
wie zum Beispiel Muskelfasern, kunstlich erzeugten Fasern und hieraus gewonnenem Gewehe sowie von kunstlichen 65 
oder naturlichen porosen Medien (worunter zum Beispiel auch Gcfafie und KapiUare in tierischen und pflanzlichen Or- 
ganismen verstanden werden sollen). 
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Paten tanspruche 

1 . Bildgebungsvcrfahren, bei dem aus Messdaten wenigstens ein Biki erzeugt wird, dadurch gekennzeichnet, 
dass Koordinaten wenigstens eines Punktes erniittelt werdcn, fur den die Konzentration von wenigstens einem Stoff 

5 einen extremaien Wert aufweist, 

dass eine zeilliche Vcriinderung der Konzentration fur mehrere Etemente eines ersten Bildes und eines weiteren Bil- 
des erniittelt wird, wobei eine zeilliche Veranderung der Konzentration C des Stoffes im Wesenllichen gem&B der 
Forrnci 

l0 dc/dt - v - (evc) 

erniittelt wird, wobei V den Nablaoperator (in 3D: (d/dx, didy, 3/dz) T ) bczeichnei und wobei D einen ortsabhangi- 
gen Diffusionstensor bezeicbnet. 

2. Verfabren nach Anspruch 1, dadurch gekennzeichnet, dass die Gleichung 
dC/dt = V • (QVC) 

fur verschiedene Zeiten gelds t wird, wobei fur wenigstens eine erste Zeit t 0 eine Ausgatigskonzenlration beruck- 
sichtigl wird, die ein ausgepragles Maximum aufweist, 
20 Die Ver.teikingsfunkl.ion wird dabei in der Umgcbung des Maximums zuihindest nafterungsweise als radialsymmc- 

triscb betrachtet. Als besonders zweckmafiig erweist es sich, wenn die Verteilungsfunktion im Wesenllichen der 
Forme I 

Co - exp(~((x - x 0 ) 2 + (y - y 0 ) 2 + (z zo) 2 )/(2a 2 )) 

entspricht, in der G die Breite der Vertetlungsfunktion angibt, die zweckmaBigerweise kleiner oder gleicb der Lan- 
25 gendimension eines Voxels gewahit wird und bei der x 0 , y 0 und Zq die karicsischen Koordinaten des (Peak-) Maxi- 

mums angeben, das an den Punkt, fur den die Konneklivitat untersucht werdcn soil, oder an einen benachbarten 
Punkt. gelegt wird. 

3. Verfahren nach einem oder beiden der AnsprCiche 1 oder 2, dadurch gekennzeichnet, dass die Komponenten der 
DiffusionsLensoren, die sich in Bereichen mil nahezu freier Diffusion, also grofien Werten der Spur des Tensors be- 

30 finden, auf einen sehr kleinen Wert (z. B, null) gesetzl werden. 

4. Verfahren nach einem oder mehreren der vorangegangenen Anspruche, dadurch gekennzeichnet, dass die Kom- 
ponenten der Diffusionslensoren, die sich in Bereichen nicdriger Diffusion- An isotropic befinden, auf einen sehr 
kteinen Wert (z. B. null) gesetzt werden. 

5. Verfahren nach einem oder mehreren der vorangegangenen AnsprCiche, dadurch gekennzeichnel, dass die Kom- 
35 ponenlen der Diffusionstensoren mil. einem (skabren) Diffusions-Anisolropie-Index, vorzugsweise mil. dem RA-In- 

dex gewichtet werdcn, wobei eine vorhcrigc Normierung der Diffusion s-Tensoren zweckmafiig is*, die vorzugs- 
weise dadurch erfolgl, dass die Elemenle eines Tensors dutch dessen Spur dividiert werden. 

6. Verfahren nach einem oder mehreren der vorangegangenen Anspruche, dadurch gekennzeichnet, dass Linien 
oder Flachen konstanter Konzentration (Isolinien, -Aachen) bestimmlund Abweichungen von einer Kreis- bzw. Ku- 

40 gel form erniittelt werden . 

7. Verfahren nach einem oder mehreren der vorangegangenen Anspriiche, dadurch gekennzeichnel, dass in jedem 
Zeitschritt die Konzentrationswerte der isolinien bzw. -flachen so gewahlt werden, dass sie einem konsianten 
Bruchleil des momentanen Maximalwertes der Konzentrationsftmktion oder des aktuellen Wertes der Konzentrati- 
onsfunktion am Ort des anfanglichen Konzenlralionsmaxiniums entspricht. 

45 8. Verfahren nach einem oder mehreren der vorangegangenen Anspriiche, dadurch gekennzeichnet, dass die Isoli- 

nien bzw. -flachen so festgelegt werden, dass sie Punkte mitderselben Krummung der Konzenlrationsfunktion ver- 
binden, (und nicht Punkte gleicher Konzentration), wobei die Krummung im Wesenllichen definiert ist als k c = 
(grad(C)) T H(C)(grad(C)), wobei H die Hesse-Matrix zur Konzentrationsfunkdon C bezeichnet (11^ - 5 2 C5x^Xj) 
und im zweidimensionaten Fail speziell gilt; 

so a 0 ' r = {dadxi, dc/B^Meaax^ + (ac/2x 2 ) 2 y 1/2 

h= {(d 2 c/d^ f d 2 c/d^dx 2 ) t (Pc/a^dx*, a 2 c/dx 2 2 )) 

-> T -> 

55 

k c = a 0 H a 0/ 

9. Verfahren nach einem oder mehreren der vorangegangenen Anspruche, dadurch gekennzeichnet, dass Punkte 
der Trajektorien dadurch erniittelt werden, dass auf den Isolinien/- flachen fur aufeinanderfolgende Zeilpunkte je- 

60 weils der Punkt ausgewahlt wird, der den maximalen Abstand zutn Ursprung aufweist, wobei zwischen den so er- 

mittelten Punkten mitlels steliger Funktionen, vorzugsweise Geraden, interpoliert wird. 

10. Verfahren nach einem oder mehreren der vorangegangenen Anspruche, dadurch gekennzeichnet, dass die Be- 
rechnung der Trajektorien sliickweisc erfolgt, wobei nach jedem Teilstuck mit einer neuen Ausbreitungsberechnung 
begonnen wird und zwei benachbarte Berechnungsvolumen einen nicht- verschwindenden tJberlapp aufweisen. 

65 1 1 . Verfahren nach einem oderen mehreren der vorangegangenen Anspruche, dadurch gekennzeichnel, dass fur die 

IsoUnicn oder Isoflachen Krlimmungen berechnet werden, wobei fiir kartesische Koordinaten die Krummung von 
Lsolinien, die durch Funktionen y(x) beschrieben werden, im Wesenllichen gemaB der Formcl 
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d 2 y / dx 2 
k = 

[ 1 + (dy/dx) 2 3 3/2 

berechnet wird. Im Pall von Isoflachen, die durch Funktionen z(x, y) beschrieben werden, erfolgt fur die Punkte (x, 
y) die Bcrechnung zweier Krummungsradien, wozu eine quadrailsche Gleichung gelost wird, die im Wesentlichen 
die folgende Form hat: 

(it - s 2 )R 2 + h[2pqs - ( 1 + p 2 )t - (1 + q 2 )r)R + h 4 - 0 
mil 

p = dtjdx, 

q - 321 By, 

s = dhjdxdy, 
t = dhJdy 1 
und wobei 

h = /i -§- p- + q 2 und aus diescn Krummungsradien eine Bcrechnung einer mitileren Krummung, im WeseniHchen 
gcmaB der Formel 

K = V4(l/R| +1/R 2 ) 
erfolgt, 

12. Verfahren nach einem oder mehreren der vorangegangenen Anspruche, dadurch gekennzeichnet, dass von fur 
aufeinander folgende Zei I punkte berechnete Isolinien oder isoflachen jene ermitteit. werden, bei denen erstmals po- 
sitive Werte fur die Krummungen der Isolinien (k) oder Isoflachen (K) auftreten (allgemciner: Wechsel des Vorzei- 
chens gegenuber dem Krtimmungswert der urspriinglichen Konzentrationsvcrteiiung nach Subtraktion einer Kon- 
stanien). Zur Vermeidung von Rauschartefaktcn kann es sich ais xweckmaSig erweisen 7 zu fordern, dass ein Tnte- 
gralwert, der Krunimungsfunktion einen bestimmten Schwcllwert uberschreitet. 

13. Verfahren nach einem oder mehreren der vorangegangenen Anspruche, dadurch gekennzeichnet, dass fur min- 
destens einen Punkt mindestens eincs Rasters mindestens ein mitllcrer Abstand mindcstens einer Isohnie oder Tso- 
fiachc vom Oil des ursprunglichen Konzentrationsmaximunis ermitteit wird, bei dem (erstmals) ein Vorzeichen- 
wechsel in k b/.w. K auflriu. 

14. Verfahren nach einem oder mehreren der vorangegangenen Anspruche, dadurch gekennzeichnel., dass Punkte 
von Isolinien oder Isoflachen, bei denen erstmals eine positive Kriimmung k bzw. K aufgetreten ist, oder benach- 
barte Punkte als Staripunkte neuer Aushreitungsberechnungen gewahlt werden. 

15. Verfahren nach einem oder mehreren der vorangegangenen Anspriiche, dadurch gekennzeichnet, dass aufein- 
anderfolgende Startpunkte der Aushreitungsberechnungen mileinander durch stetige Funktionen, vorzugsweise Ge- 
radcn, vcrhunden werden. 

16. Verfahren nach einem oder mehreren der vorangegangenen Anspruche, dadurch gekennzeichnet, dass die Fest- 
legung mindestens eines Anfangspunkts (mit den Koordinaten x 0 , yo und z 0 ) so erfolgt, dass sich dieser auf einer 
miUcls eines skaiaren, d. h. eines auf skalare Diffusions- Anisotropic- Daten angewcndeten Tracking verfahrcns, be- 
stimmten Trajektone befindet. Als Anisotropic- Daten dienen vorzugsweise fur den zugrundeliegenden Diffusions- 
tensor-Dai en satz berechnete Werte <les Anisofropie~Indcx RA. 

17. Verfahren nach einem oder mehreren der vorangegangenen Anspruche, dadurch gekennzeichnel, dass das ska- 
lare Tracking verfahren so konzipiert ist, dass es das Pfadintegral fur eine Trajektone zwischen zwei Punktcn niini- 
miert, wobei iiber das Inverse eines Anisotropie-Indexes (vorzugsweise den RA-Index) intcgrierl. wird. 

18. Verfahren nach einem oder mehreren der vorangegangenen Anspruche, dadurch gekennzeichnel, dass Zeiien 
ermitteit werden, die angeben, wann in verschiedenen Punkten bestimmte, relative Konzentrationswerte iiberschrit- 
ten worden sind und dass <iiese Daten in Time-of-ArrivaUMaps zusammengefasst werden. 

19. Verfahren nach einem oder mehreren der vorangegangenen Anspruche, dadurch gekennzeichnel, dass in Time- 
of- Arrival-Maps Zeiten t anivai durch effeklive Zeit [arrival r m -i va i = i/Wivd erselzt werden, 

20. Verfahren nach einem oder mehreren der vorangegangenen Anspruche, dadurch gekennzeichnel, dass - ausge- 
hend von einem St art punkt fur eine Reihe von aufeinanderfolgenden Punkien Gradienten der 'lime- of- Arrival- 
Maps berechnet werden, wobei der jeweils nachste Punkt. als ein Punkt jener Geraden gewahll wird, dessen Auf- 
punktdem aktueilen Punkt entspricht und dessen Richlungsvektor durch den zum aktuellen Punki gehorigen Gra- 
dienten gegeben ist. 
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